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charged bosons is considered. It is shown that the screened potential is drastically dif- 
ferent from the usual Debye one. The polarization operator of photons in plasma acquires 
infrared singular terms at small photon momentum and the screened potential drops down 
as a power of distance and even has an oscillating behavior, similar to the Friedel oscillations 
in plasma with degenerate fermions. The magnetic properties of the cosmological plasma 
with condensed W-bosons are also discussed. It is shown that W-bosons condense in the 
ferromagnetic state. It could lead to spontaneous magnetization of the primeval plasma. 
The created magnetic fields may seed galactic and intergalactic magnetic fields observed in 
the present-day universe. 
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£j ! 1 Introduction 

CO 

The potential created by an electric charge in plasma is usually described by the well known 
Debye formula, see e.g. refs. [TJ [2]. The long range Coulomb potential transforms into the 
exponentially decreasing Yukawa one: 

ON ■ because the time-time component of the photon propagator acquires a constant term which 



does not vanish when the three dimensional photon momentum goes to zero. Instead of the 
vacuum fc 2 -term the inverse propagator becomes: 

k 2 ^k 2 +U 00 (k) = k 2 +m 2 D , (2) 

where e.g. for relativistic fermions [3j: 

m 2 D = e 2 (T 2 /3 + v 2 /Tr 2 ). (3) 

These results are true if the fermions in plasma are not strongly degenerate and if the charged 
bosons do not condense. The modification of the Debye screening in the case of degenerate 
fermions was studied half a century ago [3], but the impact of Bose-Einstein condensate (BEC) 
on the screening of impurities was considered only very recently [6]. Surprisingly in the 
presence of BEC the screened potential drastically changes and becomes even an oscillating 
function of distance, which decreases either exponentially or as a power of r. 
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In the first part of this talk I discuss the impact of BEC of a charged scalar field on the 
screening of charged impurities in plasma. It is based on the works done in collaboration with 
A. Lepidi and G. Piccinelli [6]. The second part about ferromagnetism of the condensate of 
charged vector bosons is based on paper [7] of the same group. More detailed list of references 
can be found in these papers. 

Let us consider electrically neutral plasma with large electric charge density of bosons which 
is compensated by charged fermions. Bosons condense when their chemical potential reaches 
the maximum allowed value: 

HB =m B - (4) 

This can be easily seen from the consideration of the kinetic equation. Indeed the equilibrium 
distribution of bosons, if and only if their chemical potential is equal to their mass, takes the 
form: 

= CS® (q) + V , (5) 

exp l(E - mB)/T\ - 1 

where the constant C is the amplitude of the condensate. One can check that fg 9 ^ annihilates 
the collision integral for an arbitrary C. Thus, by definition, it is an equilibrium distribution. 
It is worth noting that the equilibrium distributions are always functions of two parameters: 
temperature, T, and chemical potential, fi, when fj, < m, and temperature T, and the amplitude 
of the condensate, C, when [i = m. 

We calculate the time-time component of the photon polarization operator in a simple 
straightforward way perturbatively solving operator equation of motion for the electromagnetic 
field (Maxwell equations) and averaging them over medium. One can reach the goal without 
applying to more refined real or imaginary time methods (for a review of these methods see e.g. 
Ref. [2]). 

2 Maxwell equations in thermal bath 

The standard Lagrangian of QED with interacting electromagnetic field and charged scalar and 
fermion fields with masses ms and rap respectively and with opposite electric charges ±e has 
the form: 

C = -\f^F» u - m\\^\ 2 + Kfy + i eA^\ 2 + -eJfi.- m F )^. (6) 

The equations of motion for the electromagnetic and charged scalar and spinor fields are 
respectively: 

(i$ — m)ip(x) = ejfiip(x) , (7) 

(d^ + m 2 )^) = J^x), (8) 

d u F^{x) = J»(x), (9) 
where the currents J can be written as: 



J^ix) = -ie 



d^{x) + 2A^x)d» 



{x) + e 2 A^{x)A ll (x)(l)(x), (10) 



jv{x) = -ie 



{4>\x)d^(t>{x)) - (d»^(x))(f)(x) 



+ 2e 2 A^{x)\(t){x)\ 2 - eip(x)^ip(x). (11) 
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Here F^ u = d fJi A u — d u A^ and (fTTj) is the total electromagnetic current of bosons and 
fermions. 

Operator equations (|7|) and (jHJ) can be formally solved as: 



<j>{x) 



</>o(x) + J d 4 yG B (x-y)J^y), 
%f) (x) + / d 4 yG F (x - y)e4(y)ip(y) , 



(12) 
(13) 



where Gb and Gf are the Green's functions of bosons and ferminos respectively and the zeroth 
order fields satisfy the free equations of motion: 



(8^ + m|)0 o (x) = 0, (i# - m F )Mx) = 



and are quantized in the usual way: 

d 3 q 



ip (x) 



d 3 q 



a(q)e- iqx +6 t (q)e 



iqx 



v/(2^V E 



m F 



c r (q)u r (q)e iqx + 4(qK(q)e 



iqx 



(14) 

(15) 
(16) 



Substituting solutions (|12ll3p into eq. ([9]) we obtain the Maxwell equations with the lowest 
order corrections to the electromagnetic current: 



duF^(x) 



-i e 



d 4 yG B (x - y)Jfo{y) 
d 4 yG B (x - y)J^ (y) 



eipoix^ipoix) 



- ie(j)l{x)d^ 
+ ied^(j)l{x) 

- e^ (x)j fl I d A yG F (x - y)e4(y)ip(y) 



ie 



d 4 yG B (x - y)J<t, Q (y) 



it 



it 



d yG B (x- y)J<p (y) 



d A yMy)4(y)G* F (x-y) 



+ 2e 2 A»(x)\Mx)f 



i^Mx) 

(17) 



To derive the Maxwell equations with the account of the impact of medium on the photon 
propagator we have to average operators <fi and ip over the medium. The first term in eq. 
()17p . linear in e, is non-zero if the medium is either electrically charged or possesses an electric 
current. We assume here that this is not the case, i.e. the medium is electrically neutral and 
"current-less" . 

The products of creation-annihilation operators averaged over the medium have the standard 
form: 



(at(q)a(q'))=/ B (^)5 {3) (q-q / ), 

(a(q)at(q')) = [l + /B(i? P )]<5 {3) (q-q / ), 

(c t (q)c(q'))=/F(^ (3) (q-q / ), 
( C (q)ct( q ')) = [l-/ F (i? p )]^(q-q') 



(18) 
(19) 
(20) 
(21) 



where f F)B (E q ) is the energy dependent fermion/boson distribution function, which may be 
arbitrary since we assumed only that the medium is homogeneous and isotropic. We also 
assumed, as it is usually done, that the non-diagonal matrix elements of creation-annihilation 
operators vanish on the average due to decoherence. For the vacuum case f F ^ B {E) = and 
we obtain the usual vacuum average values of (aa^) and (a^a) = 0, which from now on will be 
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neglected because we are interested only in the matter effects. As a result we obtain linear but 
non-local equation for electromagnetic field A fl (x), for which it is convenient to perform the 
Fourier transform: 



A"{k) 



d x 
(2^)3 



-ikx 



(22) 



Finally we find that field A^(k) satisfies the equation 

[k p k p g^ - k»k u + IF"(Jfe)] A u (k) = J"(k), 



(23) 

which is equivalent to the photon equation of motion (|17p in momentum space. 

In this way the photon polarization tensor, which contains contributions from the charged 
bosons and fermions, H fMU (k) = !!-,(&) + H^ u (k), according to eq. (fT7|) . can explicitly found in 
the lowest order in e 2 : 



d 3 q 



(2vr) 3 £ 



[f B (E) + f B (E)} 



1 (2g - k) ^2q - k\ 
2 



{q-kf 



m, 



l(2q + k) fl (2q + k) l 



{q + kf 



m 



B 



(24) 



2e 



d 3 c 



[f F (E) + f F (E)} 



Qy{k + q)^ - q p k p g^ v + q p (k + q) v 



+ 



(2n) 3 E 

Qu(q - k)^ + q p k p g pu + q p (q - k) u 



(k + qf 



m. 



{k-qf 



m. 



(25) 



The static properties of the medium are determined by the time-time component of the polar- 
ization tensor in the limit of uj = 0, which can be easily calculated from the above expressions: 



n ooO) - 



Hoo (k) — 2 , 

7T JO 



dqq 2 
E 

dqq 2 



C/ fl + / fl )(l + fLln 
E 2 



2q + k 



2q-k 
2q + k 



2q — k 



(26) 
(27) 



Here and in what follows k and q are respectively the absolute values of the spatial component 
of the photon and the charged particle momenta. Expressions (|26|27|) coincide with the well 
known ones found by other methods. Our new results for screening come from an addition of 
the condensate term to /b, eq. (0). 

After straightforward calculations we find that the time-time component of the charge boson 
contribution into the photon polarization tensor at zero frequency and small k (but high T) has 
the form: 



n, 



00 



[k 2 + e 2 (m 2 + m\/k + rn^/k 2 )] 



(28) 



where 



m 2 
m\ 

4 

m 2 



2T 2 /3 + C/[(2vr) 3 m B ] 

m|T/2 

4Cm B /(2vr) 3 . 



(29) 



In fact the same dependence on k is true for any temperature but the coefficients rrij may be 
different. 
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3 Screening of electric charge 



The screened potential is determined by the Fourier transformation of the photon propagator 

(^-rioo)- 1 : 

f d 3 k exp(ikr) qQ f°° dkk 2 sin At 
U(r) = Q / — ^ , 2 - ... = \ , 2 - , — — . (30) 



(27r) 3 k 2 - n o(fc) 2vr 2 7 k 2 -Il 00 (k) kr 

If IIoo is an even function of k, as is usually the case, the integration over k can be extended 
to the interval from — oo to +00 and the integral can be taken as a sum over residues of the 
integrand. For example, if the term proportional to 1/k can be neglected (low temperature 
case), Hqo is evidently even and its poles can be easily found: 



\ 



e 4 m,Q „ 4 
— — - — e z m,2 



iiem 2 ,. (31) 



The last approximate equality is formally true in the limit of vanishingly small e. 

Thus in the presence of the charged Bose condensate the "Debye" poles acquire the non-zero 
real parts: 

kj = ±^fem 2 exp(±i7r/4) = k'j + ik", (32) 
Non-zero k' leads to the oscillating behavior of the potential 

tti \ ^exp(- v ^e/2m 2 r)cos( v r e/2rn 2 r) 

U(r)j~Q . (33) 

If the term proportional to 1/k is present in IIoo, the calculations of the potential are slightly 
more complicated. Now the integration path in the complex k plane cannot be extended to 
—00 but the integration should be done along the real axis from to 00, then along infinitely 
large quarter-circle, and along the imaginary axis from 00 to 0. The result would contain the 
usual contributions from the poles and other singularities (see below) in the complex /c-plane 
and the integral over the imaginary axis. The former gives the usual exponentially decreasing 
potential, while the latter gives a power law decrease: 

U ~ Qm\/{e 2 mlr Q ). (34) 

Notice that the potential is inversely proportional to the electric charge squared. This is because 
we consider asymptotical behavior of the screened potential at large distance, r, when the 
parameter er formally tends to infinity. 

There are some other singularities in the integrand of eq. (|30p , which arise from the pinching 
of the contour of integration over q in eq. (f2l?j) by the poles of the distribution function, fs ©, 
and the branch point of the logarithm. The induced singularities in the complex k-plane give rise 
to the screening effects analogous to the Friedel oscillations [3]. If the first "pinch" dominates, 
the screened potential is: 

where z = 2r^2irTmB- Notice that U±(r) is inversely proportional to e 2 and formally vanishes 
at T — > 0, but remains finite if \JTwiBr 7^ . 

If all pinches are comparable, the screened potential drops down as a power of distance: 

U ^ ~ ~ 2e 2 T 2 m|r 6 In 3 (V8m B Tr) ' ^ 
More details can be found in the second paper of Ref. [6]. 
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4 Condensation of vector bosons 



The condensation of the charged vector bosons [7] of the electroweak group might take place in 
the early universe if the cosmological lepton asymmetry was sufficiently high. Condensation of 
vector particles differs from that of the scalars due to the additional degree of freedom, their 
spin states. Depending upon the interactions between the spins, they can be either aligned or 
anti-aligned. These states are called respectively ferromagnetic and anti-ferromagnetic ones, 
see e.g. Ref. [9j. We show that VF-bosons of the minimal electroweak theory condense in 
the ferromagnetic state and spontaneous magnetization of the primeval plasma could generate 
strong primordial magnetic fields on macroscopically large scales. 

Recently somewhat similar problem of condensation of deuterium nuclei in astrophysics has 
been studied in Ref. [8]. The authors argue that the interaction between deuterium nuclei forces 
them into the lowest spin antiferromagnetic state. 

In the minimal standard electroweak model the spin-spin interaction of VF-bosons is deter- 
mined by the interaction between their magnetic moments and their contact quartic coupling. 
The former can be found from the analogue of the Breit equation for vector particles which 
leads to the spin-spin potential of the form [7J: 



UZ n (r) = ^ - 3 (Sl ' r) 5 (S2 - r) - *%- ( Sl • S 2 ) ^)(r 



rn 



(37) 



where a = e 2 /4ir ~ 1/137 and p is the ratio of the real magnetic moment of W to its value in 
the electroweak theory. Since the plasma is supposed to be neutral, the Coulomb interaction 
between the condensed W is compensated by the charged leptons. 

To find the energy shift of a pair of VF-bosons due to this interaction we need to average 
potential (|37p over the W-wave function. The wave function of the condensate is supposed to 
be angle independent S-wave state. Thus the energy shift induced by the magnetic spin-spin 
interaction is expressed through the integral of potential (f37|) over space: 

/£<«•('> --^(fc-s,), < 38 > 

where V is the normalization volume. 

Since S 2 ot = (Si + S2) 2 = 4 + 2S1S2, the average value of S1S2 is equal to 

S 1 S 2 = S 2 tot /2-2. (39) 

For Stat = 2 this term is S1S2 = 1 > 0, while for Stot = it is S1S2 = —2 < 0. Thus, if 
the spin-spin interaction is dominated by the interactions between the magnetic moments of W 
bosons, the state with their maximum total spin is energetically more favorable and VF-bosons 
should condense in the ferromagnetic state. Such an interaction would lead to the spontaneous 
magnetization of the vector particles in the early universe. 

Another contribution to the spin-spin interactions of W comes from their quartic self- 
coupling: 



e 2 



2 



r 1 e 2 fw t x W^. 

L «" = [K^w„w» - (wtw'f] = \ sinigw . (40) 



It is assumed here that d^W^ = and thus only the spatial 3- vector W is non- vanishing, while 
Wt = 0. The Fourier transform of this term with proper (nonrelativistic) normalization leads 
to 

^ = 2 e2 2 fl (SiS 2 )^(r). (41) 
8m|^ sm z 9 W 
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Thus the quartic self-coupling of W contributes only to the spin-spin interaction whose sign is 
antiferromagnetic. 

In the minimal standard model the interaction between the magnetic moments of W domi- 
nates and the condensed VF-bosons have ferromagnetic behavior. However, in some modification 
of the standard model antiferromagnetic behavior is possible and W would condense in the state 
with zero or microscopically small total spin. In this case classical vector field of VF-bosons could 
not be created, though they would still make the Bose condensed state. 

The exchange of Z°-bosons may also contribute to spin-spin interactions of W. It can be 
shown that for non-relativistic Z this contribution vanishes [7]. However, if the momentum 
carried by the virtual Z is non-negligible in comparison with its mass, the contributions of Z 
and photon exchanges are similar. 

The long range interactions between the magnetic moments of TV-bosons, in principle, can 
be screened by the plasma effects. However, in contrast to electric interactions, which are Debye- 
like screened, magnetic interactions in pure electrodynamics (or in any other Abelian theory) 
are known to remain unscreened. On the other hand, in non-Abelian theories the screening 
may occur in higher orders of perturbation theory due to the violent infrared singularities, 
which make impossible perturbative calculations [10J . The screening may potentially change 
the relative strength of the electromagnetic spin-spin coupling, which is affected by screening 
effects, with respect to the local quartic, iy 4 -coupling which is not screened. However, in the 
broken phase the system is reduced to the usual electrodynamics, where screening is absent and 
IV-bosons would condense in the ferromagnetic state. In the unbroken phase of the electroweak 
theory the answer is not yet known. 

If the magnetic intraction between the spins of W-bosons is screened the potential describ- 
ing the magnetic spin-spin interaction is related to amplitude (|37p with a modified photon 
propagator. So it can be written as: 

Ut m) (r) =-^[ ( % 1 Tl^L W (Si • S 2 ) - (Si • q) (S 2 • q)l , (42) 

m w J ( 27r ) W + I Mq)) L J 

where H ss is the plasma correction to the space-space component of the photon propagator. 

If, as above, we assume that the wave function of IV-bosons is space independent and average 
this potential over space, we obtain the following expression for the spin-spin part of the energy 
shift: 

" - / P^iO -^rj f^^LM m 

Clearly SE vanishes if U ss is non-zero at q = 0. Of course, this is an unphysical conclusion, 
because the integration over r should be done with some finite upper limit, r max = I, presumably 
equal to the average distance between the W bosons. So instead of the delta- function, 5^ 3 ^(q), 
we obtain: 

rl 4^ 

/ d 3 r exp(iqr) = [sm(ql) — qlcos(ql)] . (44) 
Jo <T 

The energy shift is given by the expression: 

r P _ A eV gigj f d 3 q [sin(gl) - ql cos(gl)] [g 2 <% - gV] 

5E ~ ~^v^ s ^ 1 Wf TWTTiZM] ' (45) 

where V = 4vr/ 3 /3. 

When we average over an angle independent wave function, e.g. S-wave for the condensate, 
the non-vanishing part of the integral in Eq. ()45j) is proportional to the Kronecker delta, hence: 

5E = S[S 3 2 A5 ij , (46) 
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where the coefficient A can be calculated by taking trace of Eq. f|45|) : 



Hence the energy shift of a pair of VF-bosons in S-wave state due to the spin-spin interaction 
is: 

- f a a ^ 87re V /" A [ sin (gQ ~ cos(gZ)] 

^ " " (Sl • S2) 3^4 7 (2^3 q tf + n s M] • (48) 

Introducing the new integration variable x = ql, we can rewrite it as: 

^ = - (Sl ' S2) 8^4/0 ^TiXp) [^mx + / 2 n s ,(x//)cosxj , (49) 

We used here the usual regularization of divergent integrals: exp(±iql) — > exp(±iql — eq) with 
e — > 0. With such regularization Jq 00 <ixcos(x) = 0. 

Evidently, if II SS = 0, we obtain the same result as that found above. In fact the necessary 
condition for obtaining the "unscreened" result is l 2 TL ss (x/l) <C 1, but for a large Z 2 II SS the 
electromagnetic part of the spin-spin interaction can be suppressed enough to change the ferro- 
magnetic behavior into the antiferromagnetic one. This might take place at high temperatures 
above the EW phase transition when the Higgs condensate is destroyed and the masses of W 
and Z appear only as a result of temperature and density corrections and thus are relatively 
small. The quantitative statement depends upon the modification of the space-space part of 
the photon propagator in presence of the Bose condensate of charged W. As far as we know, 
this modification has not yet been found. 

If VF-bosons make a ferromagnetic state, the primeval plasma, where such bosons condensed 
(possibly due to a large cosmological lepton asymmetry), could be spontaneously magnetized, 
as it happens in usual ferromagnets. The typical size of the magnetic domains is determined by 
the cosmological horizon at the moment of the condensate evaporation. The latter takes place 
when the neutrino chemical potential, which scales as temperature in the course of cosmological 
cooling down, becomes smaller than the W mass at this temperature. 

The large scale magnetic field, produced by the ferromagnetism of W^-bosons, might survive 
after the decay of the condensate due to the conservation of the magnetic flux in the primeval 
plasma because of its high electric conductivity. Such magnetic fields, which were uniform 
at macroscopically large scales at formation, may be the seeds for the observed galactic or 
intergalactic magnetic fields at astronomically large scales. Evidently the characteristic size of 
magnetic domains at production is much smaller than the galactic size, even with an account 
of the cosmological stretching out. Nevertheless, magnetic fields homogeneous at astronomical 
scales may be created by chaotic reconnection (Brownian motion) of the magnetic field lines 
at much smaller scales but by an expense of the field amplitude. Such a mechanism could be 
a competing alternative among many other attempts to solve the mystery of the generation of 
large scale magnetic fields in cosmology, for a review see ref. 
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